typeset using JPSJ.sty <ver.l.Ob> 



Numerical Study on Aging Dynamics in the 3D Ising Spin-Glass Model. 
I. Energy Relaxation and Domain Coarsening 
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Using Monte Carlo simulations, we have studied the relaxation of energy of the three- 
dimensional Ising spin-glass model in aging process. Our finite-size-scaling analysis on the 
isothermal energy decay after the quench suggests strongly that the energy relaxes by coarsening 
of domain walls in agreement with the droplet theory. We have also evaluated relaxation times 
required for spin configurations in small systems, which are regarded as isolated droplets, to flip 
globally. The results support the fundamental assumption of the droplet theory that coarsening 
of domain walls, in case of spin glasses, is driven by successive nucleation of thermally activated 
droplets. 
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§1. Introduction 

In recent years ofF-equilibrium dynamics, particularly 
aging dynamics, of stdn glasses Iiels attracted much 
attentioiL theoretically,EP numericallyp' and experimen- 
tally.EroEl' Many experiments have confirmed that the re- 
sponse to small dc/ac magnetic field after quenching the 
system below the spin-glass transition temperature 
shows aging effects persistently at least up to the largest 
possible time scale available in laboratories. On the side 
of theories, oneLoLthe most notable advances is the an- 
alytical theorjoBLO-* based on solvable mean-field spin- 
glass models, which has provided non-trivial predictions 
on aging effects. 

Presumably one expects that the slow relaxation is 
caused by some underlying collective spin excitations in 
the system. In the case of ordering processes in non- 
random systems like the conventional Ising ferromagnet, 
the relaxational dynamics after quenching the system be- 
low the phase transition temperature is related to coars- 
ening of domain walls which separate ordered domains 
of two different pure states .cP The typical size of the do- 
mains, R{t), at time t usually grows by a certain power 
law. This behavior can be directly studied, for instance 
by scattering experiments. The intriguing question is 
if such a picture of the coarsening of domain walls also 
holds for aging processes in spin glasses. Obviously, the 
above mentioned mean-field theory, which is exact only 
in infinitely large dimension, is not helpful in this respect. 
TVip dnrnajn picture based on the so-called droplet the- 
ory,Q'E3EJ'Ej'' has provided an interesting phenomenol- 
ogy. However the pure states of spin glasses cannot be 
found in practice which makes it difficult to rationalize 
the domain picture in a-jdirect way. 

In the present studyji3) we mainly focus on the energy 
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relaxation in the 3-dimensional EA spin-glass model dur- 
ing the aging process simulated by Monte Carlo dynam- 
ics. The latter starts from arbitrary initial configura- 
tions at temperature T below which simulates instan- 
taneous quench to T from T = oo. We have found that 
the relaxation of energy can be explained consistently in 
terms of coarsening of domain walls as in the case of usual 
phase ordering processes in non-random systems. Let us 
briefly describe the idea here. Suppose a system is aged 
by elapsed time t after the quench, and a typical size of 
domain walls has grown to R{t). In the presence of do- 
main walls, the energy per spin, exit), is larger than the 
equilibrium energy e^'* at temperature T by an amount 

deT[t) = eT[t) - e>p ' (X — . , . , , , — . (1.1) 



with ^0 being a certain microscopic unit of length. Here 
the numerator on the right-hand side is the domain wall 
energy where T(T) is the stiffness constant and 9 is the 
characteristic exponent. The denominator is the volume 
of a domain with d being the dimension of the space. 

We have studied systematically the energy per spin on 
relatively small system sizes L and at relatively higher 
temperatures in the spin-glass phase. This enables us to 
see the behaviors of two time regimes within our compu- 
tational time window; the size-independent behavior of 
domain growth in the dynamic regime R{t) < L and the 
size-dependent equilibrium behavior in the static regime 
L < R{t). It has turned out that the simulated data in- 
cluding in the crossover region between the two regimes 
can be concisely described by a finite-size-scaling (FSS) 
function. By multiplying eq. (1.1) by a scaling function 
Se{R{t)/L), the scaling ansatz is proposed as 



Serit) = T(r) (^) 



(1.2) 



with R{t) = ti/^(^). The FSS analysis allows us to obtain 
the characteristic exponent 9 for the domain wall energy. 
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The value of 6 turns out to be nearly equal to the result 
by Bray and Moore at T = oB-' These results strongly 
support that the aging process in the present spin-glass 
model is described well in terms of coarsening of domain 
walls, i.e., the droplet picture. 

As described just above the FSS analysis on the en- 
ergy relaxation yields the power-law growth in t of the 
typical domain size R{t) with the exponent z{T), which 
is consistent with that obtained tg£_thn4jrevious numer- 
ical studies using replica-overlapJl3' oils' We have also 
performed analysis on the replica-overlap with greater 
statistical accuracy, and have established the consistency 
between the two analyses. 

A fundamental assumption of the droplet theory on 
aging processealj-' is that coarsening of domain walls, 
in case of spin glasses, is driven by successive nucle- 
ation of thermally activated droplets. In order to test 
this assumption, we have studied largest relaxation times 
tl{T), which are those needed for a global spin- flip in 
each sample with relatively small sizes L. We have found 
that the L-dependence of tl{T) is consistent with the 
growth law of domains, i.e., R{t ^ tl{T)) ~ L. Fur- 
thermore, the temperature dependence follows in fact an 
Arrhenius law. 

In our simulated results, however, there remain some 
discrepancies with the original droplet theory. Firstly, it 
is assumed to grow as L"^ in the droplet theory, while 
the barrier free energy simulated grows logarithmically 
with L. This results agree with the prevkais-ponjectures 
extraxted from the i-dependencc of i?(t).Ot3' Secondly, 
the width of the distribution of the barrier free energy 
does not increase with L. 

The present paper is organized as follows. After de- 
scribing the model and numerical method in the next 
section, we present and discuss the results of our sim- 
ulation in § ^ Section ^ is devoted to the concluding 
remarks on the present work. 

§2. Model and Numerical Method 

We carry on simulations on the 3D EA Ising spin-glass 
model with nearest-neighbor interactions {Jij} defined 
by the Hamiltonian, 



H 



JijSiSj, 

iv) 



(2.1) 



on an N = L X L X L simple cubic lattice. The periodic 
boundary conditions are adopted. The quenched random 
interaction {Jij} are drawn from a Gaussian distribution 
with zero mean and variance one. It has been reported 
that this model has a spin-glass ordered phase below the 
critical temperature Tc whose most recent value is re- 
ported as Tc = 0.95 ± 0.04.0 We use the single-spin-flip 
heat-bath Monte Carlo method. As a unit time, which 
we call one Monte Carlo Step (MCS) hereafter, we adopt 
the time in which N spins are updated. 

The energy per spin exit) at t MCS after the temper- 
ature quench is calculated as 



1 

TV 



1 



t + tta 



2tt 



= t-tta 



H( 



(2.2) 



where H{t) is a value of Hamiltonian at r MCS. In order 
to reduce thermal fluctuations we take here an average 
over short time 2tta + 1 around t with tm — t/lOOO. The 
bracket [• ■ -Jav indicates the average over Ng Monte Carlo 
runs with independent realizations of quenched random 
variables {Jij}, initial spin conflgurations, and random 
numbers. We have simulated systems with relatively 
small sizes (L = 4 ~ 10) in order to study the crossover 
between the dynamic and static regimes as mentioned in 
§1, as well as large systems with i = 24 and 32 in order 
to obtain bulk behavior which is not affected by finite 
size effects within our time windDw-/^ 10^ MCS). 

Following the previous works JIj'EJ'ES the length scale 
of ordered domains R{t) is estimated by analyzing the 
replica-overlap function, G{r,t), defined as 



G(r,t) = 



1 ^ 

i5:[5(")(t)5f)(t)5(:)(t)5if:.(t) 



(2.3) 

Here a and /3 are indices of the two replicas which have 
different random spin configurations at t — (an instant 
of the temperature quench), and are updated indepen- 
dently. Its correlation length is considered to be propor- 
tional to R{t). In eq. ( |2.3| ) r is a spatial distance which 
we take only along the directions of the lattice axes. We 
have simulated large systems with L = 24 and 32 in this 
analysis. 

§3. Results and Discussions 

3.1 Finite- size scaling of exit) 

Let us begin with analysis on the relaxation of the 
energy per spin exit) of large systems with L = 32 
and 24 after the quench. The log- log plots of 6eT{t) = 

(t) - e^°°^ vs. t are shown in Fig. 1, where e^'' de- 
notes limf^ooliniip— |ooeT(0- We fitted the data to the 
following formulajla' 



(3.1) 



with c and A being the fitting parameters. From 



eq. (ITt]) with R{t) cx ^^/^(T) obtain 
A {d-e)/z. 



(3.2) 



As seen in Fig. 1, the fitting works very well. The values 
of exponent A depend on temperature. It is proportional 
to T ai^low temperatures in agreement with the previous 
workES' Below we determine z{T) and check eq. (3.2) 
numerically. 



In Fig. 2 we show the log- log plot of Serit) vs. t at 
T = 0.7 for small systems with L = 4, 5, 6, 7. It is seen 
that there exist two time regions of the relaxation curves. 
In the shorter time region the curves have no size depen- 
dence (the dynamic regime), while in the longer time 
region the curves depend on L and exit) approaches to 
an equilibrium value ei^^ which also depends on L (the 
static regime) . The crossover between the two regimes is 
expected to occur when a typical size of ordered domain 
R{t) reaches L. We therefore have tried some finite-size- 



scaling (FSS) analyses based on eq. (1.2). It has turned 
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the following asymptotic limits: 



T{T)R{tf- 
V{T)L^-'^ 



^(oo) _ J l{i )ti{t)" for dynamic regime, 

for static regime, 

(3.5) 

with R{t) being proportional to t^/'^'-'^K 

In practice, we have performed the above FSS fit in the 
following way. By making use of ei^"^ evaluated above 
using the bulk data, and expecting that Sexit) varies 
as a function of a; = \ogt and y = logL, we fitted the 
data of Sexit) to eq. (3.3) by means of the following trial 
function with ten parameters, 



\og{eT{t) 



,(oo)\ 



(0-d)logL + log/(t/L^(^))(3.6) 



Z — ^ 77 I 



Fig. 1. The power law decay of excess energy per spin (5ey(t) 
in the dynamic regime. The con stan ts are obtained by 

fitting the simulated date to eq. ( |?.]| ). The data are at T = 
0.4, 0.5, 0.6, 0.7, 0.8 from top to bottom. The systems simulated 
are with L = 32 and TVs = 1600 for T = 0.8 and 0.7, and with 
L = 24 and A^s = 3200 for T = 0.6 ~ 0.4. They are sufficiently 
large to obtain size-independent (bulk) energy decay. 




0.01 



Here we have expanded logarithm of / with respect to 
logarithm of its argument up to the seventh order. The 
above fitting has been done for 10 independent exit), 
each of which is obtained by averaging over Ns/10 sam- 
ples. We have then estimated errors of exponent and 
z{T) from its variance over the 10 independent values. 
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Fig. 2. The decay of excess energy per spin SeT{t) for systems 
with L = 4 (circles), 5 (triangles), 6 (squares), 7 (solid circles), 
and 32 (solid triangles). The asymptotic energy per spin is 
determined as —1.66204 by fitting data of L = 32 to eq. (|3.l|). 
The sample numbers A^s are 10^ for L = 4, 5, 6 , 3.5 X 10'' for 
L = 7, and 640 for L = 32. 



out that the following FSS ansatz works best: 

Serit) = exit) - 6^?°' = L'-^ /{t/ L<^^), 

^2,-A(T) Jqj, X < 1, 



T'(r) for x > 1, 



(3.3) 
(3.4) 



where T'(T) is a constant proportional to T(r) in 
eq. (1.1). In other words, the above FSS function has 



Fig. 3. The scaling plot of 5 exit) versus t/L^^'^\ 



The results of the above FSS are shown in Fig. 3. We 
see that the scaling works quite satisfactorily. The sim- 
ilar results are obtained also at T = 0.8, for which we 
have examined systems with L — 4, 6, 8, 10 and Ng up to 
130000. The exponents obtained are listed on Table I for 
and Table II for z{T). The values of ^(r) are consistent 
with those extracted from R{t) in the larger systems, as 
will be discussed in the next subsection. We therefore 
consider that the present FSS with eqs. (3.3) and (3.4) 
is quite reasonable. The exponent 9 seems almost inde- 
pendent on temperature as expected. Actually, its values 
derived by the present FSS at T = 0.7 and 0.8 are almost 
equal to the droplet exponent 6 = 0.19 ± 0.01 extracted 
from the defect energy analysis at T = by Bray and 
Moore.EP This result in the dynamic regime confirms the 
domain coarsening picture described in §1. 
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It is worth noting that static energy of finite size sys- 
tems is higher than that of the infinite one by an amount 
L^, that is, 



T'(r) ( ^ 



(3.8) 



This is interpreted that the imposed periodic boundary 
conditions distort the spin correlations withpcespect to 
what would be realized in an infinite system.113' A simi- 
lar and more transparent situation can be seen in a pure 
ferromagnet with the anti-periodic boundary condition. 
The latter induces competition or frustration in the loop 
of the interactions winding the whole system. Clearly it 
leaves a domain wall even in equilibrium. The FSS of 
eqs. (3.3) and (3.4) works well also for this case and the 
well known exponents, z — 1/2 and — d ~ 1 are de- 
duced. The situation in the present spin-glass model is 
more complicated since we cannot separate the effects of 
frustrations in bulk and of the constraint by the bound- 
ary conditions. However, the near agreement of the expo- 
nents 9 atT = 0.7 and 0.8 with the one at T = suggests 
that the periodic boundary condition alone introduces a 
'defect energy' similar to the one at T = discussed 
by Bray and Moore. In this respect let us note that a 
similar algebraic size dependence of finite size correction 
to the ground state energy with the^periodic boundary 
condition has been reported recently.113' 

3.2 Domain growth after quench 

Here we examine the replica-overlap function G{r, t) of 
eq. (2^) and the domain growth R{t) after the quench. 
Our analysis follows the work of Kisker et al£3 We use 
an improved method of evaluating R{t) from G{r, t) as 
explained below. Because of the periodic boundary con- 
ditions, relatiop^G'(r, t) = G{L — r,t) holds, and so R{t) 
is evaluated as& 



R{t)^2-—y 



G{0,t) ^ 
G{2n/L,t) 



where 



G{k,t) 



L/2 



G{r, t) cos fcrdr. 



(3.9) 



(3.10) 



If G(r, t) is simply proportional to cosh[2(r — L/2)/R{t)] 
(or exp(—2r / R{t)) in the limit L — > oo), the estimation 
of eq. (^) gives R{t) exactly. We consider that our es- 
timation of R(t) is superior to the method based on the 
moments Jq r"G{r)dr with n = 0, 1 used by Kisker et 
al. since our estimation significantly reduces the irrel- 
evant contribution from small r which deviate from an 
exponential law. 

The averaged domain size R{t) estimated from 
eqs. ( |3.9| ) and ( 3.10 ), and shown in Fig. 4, are well fit- 
ted to R{t) = bt^/^^^K The obtained exponents z{T) 
are 7.86(3), 8.71(3), 9.84(5), 11.76(6) and 14.80(7) at 
T = 0.8, 0.7, 0.6, 0.5 and 0.4, respectively. These values 
at T = 0.8, 0.7 listed on Table II are in good agree- 
ment with those extracted from the FSS analysis on 
5eT{t). The temperature dependence of z{T) is given by 
l/z(T) ~ 0.17r at T < 0.7. This linear T-dependence 
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Fig. i. The length scale R{t) of ordered domains estimated from 

eq. for temperatures T = 0.4,0.5,0.6,0.7, and 0.8 from 

bottom to top. The solid lines are fitting lines R{t) oc 




Fig. 5. Logarithms of the replica-overlap function G{r, t) versus 
at t = 5" (n = 2, . . . , 8) and r > 3 for T = 0.7 with 
L = 32 and A^s = 3500. The solid line is only a guide to the eye. 



Table I. Exponent 6 evaluated by our finite-size-scaling (FSS) 
analysis on the energy relaxation. 



T 


FSS 


Bray and Moore 


0.8 


0.20(2) 




0.7 


0.19(3) 




0.0 




0.19(1) 



of llz{T) is consistent with the resjilts by Marinari et 
al.,t3 and those by Kisker et al.,ti|) though its coeffi- 
cient 0.17 is different from the one estimated from those 
of Kisker et al. 

The replica-overlap function G{r, t) is well scaled when 
it is plotted against r/t^^'^'-'^^ as shown in Fig. 5. Fur- 
thermore it is seen that G{r, t) in a longer length scale 
decays in a simple exponential form. The exponents z{T) 
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Table II. Exponent z{T) estimated from our finite-size-scaling 
(FSS) analysis on the energy relaxation and directly from cor- 
relation length of G{r,t). The values are consisteritnwith each 
others. The values z{T) of the previous workstSEj' are also 
listed, where the values of Marinari et al. are calculated from 
their estimation as l/z(T) ~ 0.16T. The present result is con- 
sistent with z{T) by Marinari et al. but not with z{T) by Kisker 
et al. 



T 


FSS 


m 


previous works 


0.8 


7.73(7) 


7.86(3) 


~ 7.81 (Marinari ct al.) 


0.7 


8.83(8) 


8.71(3) 


~ 8.93 (Marinari et al.) 








12.4 (Kisker ct al.) 



used in this plot are those extracted by the fit of R{t). 

3.3 Distribution of barrier free-energy of isolated 
droplets 

As mentioned in §1, another important purpose of the 
present work is to study directly relaxational dynam- 
ics of droplets. Since, however, it is hard to specify 
each droplet in bulk, we examine here characteristic time 
scales tl (T) for spin configurations of systems with small 
L to turn over as a whole in equilibrium. We regard such 
systems as isolated droplets of size L. 

At first, to MCS is elapsed starting from a random spin 
configuration. Then the spin configuration {>S'|"''(io)} 
is copied to another replica /3 and the two replicas are 
updated independently with different sequence of ran- 
dom numbers. The latter is continued until the clone- 
correlation function 



1 



N 



i=l 



(a) 



{t)sl^\t). 



(3.11) 



first crosses zero at some time ti. Then the spin config- 
uration {S'^°'\ti)} is copied to rephca /3 and the above 
process is repeated until t2 when Q{t) crosses zero. The 
procedure is repeated once more, and the time when Q{t) 
crosses zero is denoted by ^3. 

In Fig. 6 we show distribution of the logarithmic times 
P£,(lnTi), where Ti — ti — with to = 0. One can see 
clearly that PL(lnr2) and ^^(Inra) coincide with each 
other, but it is not the case for P£,(lnTi). The result is 
interpreted that (In ti ) is affected by off-equilibrium 
effects, while PL(lnr2) can be regarded as an equilib- 
rium distribution which we analyze in detail below with 
Ti(T) = T2 (= r). In this analysis we have examined 
both the periodic and free boundary conditions. 

In Fig. 7-(a) we show the log-log plot of PL(lnT) versus 
r of the system with L = 4 at T = 0.8, 0.7, 0.6 and 

0. 5. If the relaxation process is of the Arrhenius-type, 

1. e., T = tq exp(i?/T), the distribution of the free-energy 
barriers ^l{B) is written as 



^l{B) - lpi(lnr/ro) 



(3.12) 



where Tq^ is the attempt frequency. With tq — 1 the 
data in Fig. 7- (a) lie top on each other as shown in Fig. 7- 
(b). This result indicates that the relaxational dynamics 




0.01 



Fig. 6. The distribution of global-flip times ti (solid circles), T2 
(solid triangles) and T3 (solid squares), at T = 0.8 in L = 4 
systems. 



of present interest is in fact an Arrhenius type whose 
barrier free energy B little depends on temperature. 

Making use of relation B = Tlnr confirmed above, 
we next examine the L-dependence of ^l{B). In Fig. 8 
we show the plot of ^l{B) against B — Bl at T = 0.8, 
where B^ is the value, at which ^l{B) is maximum for 
each L. The resultant Bl is well fitted to 



Br = A'lnL, 



(3.13) 



as shown in Fig. 9 below. Another interesting observa- 
tion in Fig. 8 is that the shape of the distribution ^l{B) 
docs not depend on L, while its peak position does. This 
is certainly in disagreement with the original droplet the- 
ory,t3' in which the width of ^l{B) is conjectured to be 
of the same order as that of its mean. 



If the coarsening process of domains represented by 
R{t) analyzed before is also assumed due to nucleation of 
domains of the corresponding size by thermally activated 
process, the associated barrier free energy Br is given 
bylif 

BR = A\nR. (3.14) 

In Fig. 9 we plot Bn thus evaluated and Bl of eq. ( 3.13| ) 
against the properly scaled length I. For Bn we set / = 
2R{t), while for B^ we use I = L/{2^/'^) since the clone- 
correlation function can cross zero when a half of spins in 
the system turns over. The logarithmic dependence on I 
of both Bft, and i?L is clearly seen in the figure. Also it is 
seen that the coefficient A for Bn is a little larger than 
A' for Bl. This difference may indicate that dynamics 
of ordered domains in a large system is not perfectly 
independent from their neighboring ones as was assumed 
in the original droplet theory. A similar conclusion has 
been also deduced experimentally.^ 
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Fig. 7. (a) The distribution of r (= rL(T)) for L = 4 at 
T = 0.8 (open circles), 0.7 (solid circles), 0.6 (open triangles) 
and 0.5 ( solid triangles) . (b) The plot of "J>i (B) extracted from 
eq. (El2). 
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Fig. 8. The plot of $i(S) versus B - B^, a.t T = 0.8 for L = 
4 (open circles), 6 (solid circles), 8 (open triangles), 10 (solid 
triangles) and 12 (open squares). 
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3.4 Discussions 

Here let us comment on the magnitude of energy, 
length and time scales in the aging process simulated 
so far. Our result is consistent with Fl ^ , while 
the barrier free energy associated with them is pro- 
portional to InL. It is then naturally expected that the 
latter is overcome by the former at a certain value of 
L, and so that the present results cannot provide us a 
proper picture in the thermodynamic limit. But we note 
that Fl evaluated from the present simulation is given 

by 



Fl oc NSct =i 2.0 x {L/lof 



(3.15) 



with 9 = 0.19, where the factor 2.0 is determined from 
the saturated value in the scaling plot of Fig. 3. Because 
of this small value of 9, combined with a relatively large 
value of A (~ A') in eq. (3.14) which in turn is due to a 
relatively large value of z{T){~ 8 ~ 10), F^ and Bn (or 
Bl) become of the same order of magnitude only when L 
is of the order of 10^ ^ 10^ with = 1- The correspond- 
ing time scale, for R{t) to reach to this length scale, is 



Fig. 9. The characteristic barrier free energy corresponding to 
characteristic length scale / (at T = 0.8). The open circles are 
estimated from the growth of R{t). The solid circles and solid 
triangles are estimated from the distribution of global flip times 
Tj^(T) with the periodic and open boundary conditions, respec- 
tively. 



much more than an astronomical one. In this context, 
it is pointed out that Bn experimentally extracted by 
Joh et alE3 recently is consistent with the simulated 
results expressed by eq. ( 3.14 ) (with t/tq ^ 10^^ and 
R/Iq ~ sevral 10). We therefore consider that the re- 
sults obtained in the present work are applicable to ag- 
ing phenomena observed within a time window of ex- 
periments on realistic spin glasses. However, mechanism 
which gives rise to the Ini? dependence of Bn is not clear 
yet. It might be related to fluctuations governed by the 
T = Tc criticality since temperatures investigated in the 
present work may not be far enough from T^. This prob- 
lem is left for future work. 
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§4. Concluding Remarks 

By means of Monte Carlo simulations we have stud- 
ied aging dynamics after the temperature quench from 
T = cxD to T in the spin-glass phase in the 3D Ising EA 
model. The results of the energy evolution erit), in par- 
ticular, the finite-size scaling on dexit) (= exit) — e^^), 
strongly suggest that the droplet picture describes ag- 
ing dynamics appropriately. It is emphasized that these 



results are obtained based on eq. (1.1) which is a most 
fundamental ingredient of the droplet picture, and that 
they are extracted by analyzing time evolution of one 
spin configuration as compared with the previous works 
which were making use of the replica-overlap. 

Also by analyzing relaxation times tl required for a 
global flip of spin configurations in small systems (iso- 
lated droplets) of size L, we have shown that the L- 
dependence of tl is consistent with the t-dependence of 
the mean domain size R{t), i.e., R{t ~ tl{T)) ~ L. This 
result further supports the droplet picture that coarsen- 
ing of domain walls is in fact driven by successive nu- 
cleation of thermally activated droplets. The explicit 
form of this t (T)-dependence of R{t) (L), however, dis- 
agrees withL.±he original droplet theory due to Fisher 
and Huse:tlHl3' the barrier free energy associated with 
ordered domains of length scale L is proportional not 
algebraically as predicted by the their theory but loga- 
rithmically to L. 
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